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ABSTRACT 

We assume that the supersymmetric hghtest neutrahno is a good candidate 
for the cold dark matter in the galaxy halo, and explore the possibility to 
produce extended diffuse radio emission from high-energy electrons arising from 
the neutralino annihilation in galaxy clusters whose intracluster medium is 
filled with a large-scale magnetic field. We show that these electrons fit the 
population of seed relativistic electrons that is postulated in many models for 
the origin of cluster radio halos. If the magnetic field has a central value of 
3 -=-30 //G (depending on the Dark Matter profile) and is radially decreasing from 
the cluster center, the population of seed relativistic electrons from neutralino 
annihilation can fit the radio halo spectra of two well studied cluster: Coma and 
1E0657-56. The shape and the frequency extension of the radio halo spectra 
are connected with the mass and physical composition of the neutralino. A 
pure-gaugino neutralino with mass > 80 GeV can reasonably fit the radio 
halo spectra of both Coma and 1E0657-56. This model provides a number of 
extra predictions that make it definitely testable. On the one hand, it agrees 
quite well with the observations that (i) the radio halo is centered on the cluster 
dynamical center, usually coincident with the center of its X-ray emission, 
(ii) the radio halo surface brightness is similar to the X-ray one, and (in) the 
monochromatic radio power at 1.4 GHz correlates strongly with the IC gas 
temperature. On the other hand, the same model predicts that radio halos 
should be present in every cluster, which is not actually observed, although 
the predicted radio halo luminosities can change by factors up to ~ 10^ 10^, 
depending on the amplitude and the structure of the intracluster magnetic 
field. Also, neutral pions arising from neutralino annihilation should give rise to 
substantial amounts of diffuse gamma-ray emission, up to energies of order M^, 
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that could be tested by the next generation gamma-ray experiments. 

Subject headings: cosmology: general - Dark Matter: neutralino - galaxies: clusters: general 
- radio emission: general - radiation mechanisms: non-thermal. 
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1. Introduction 

Dark Matter (DM) interaction and annihilation in the halo of our galaxy and in other 
galaxies have relevant astrophysical implications. In fact, if DM is constituted by weakly 
interacting massive particles (WIMP's), their annihilation can produce direct and indirect 
signals such as observable fluxes of positrons (Silk and Srednicki 1984, Rudaz and Stecker 
1988, EUis et al. 1989, Stecker and Tylka 1989, Kamionkowski and Turner 1991, Baltz 
and Edsjo 1998), antiprotons (Chardonnet et al. 1996, Bottino et al. 1998), gamma-rays 
(Bengtsson et al. 1990, Chardonnet et al. 1995) and neutrinos (Silk and Gondolo 1999) 
from the Milky Way halo. 

Motivated by these results, we explore here some specific consequences of the 
production of positrons and electrons from the WIMP interaction in massive DM halos hke 
clusters of galaxies. The decays of WIMP annihilation products (fermions, bosons, etc.) 
yield, among other particles, energetic electrons and positrons up to energies comparable to 
the WIMP mass (usually of the order of tens to hundreds GeV). It is straightforward to 
realize that these energetic electrons and positrons (hereafter we will refer to these particles 
as electrons because their distinction is not essential for our purposes) can emit synchrotron 
radiation in a DM halo which is filled with a magnetic field at the level of /iG. Clusters of 
galaxies are the largest bound systems which have indeed both the largest amount of DM 
and extended Intra Cluster (IC) magnetic field at the level of a few /xG. It is then natural to 
expect that the WIMP annihilation products can give origin to an extended radio emission. 

Extended radio halos and relics are observed, at present, in more than 40 galaxy 
clusters (see, e.g., Owen et al. 1999, Feretti 1999, Giovannini 1999, Liang et al. 2000). 
Cluster radio halos show, generally, a regular morphology, which is similar to the X-ray 
morphology, a low surface brightness, and a steep radio spectrum, J(i/) ~ i/""'', with 
~ 1 -j- 1.5. The level of polarization in the radio halo emission is usually ^ 10% (Feretti 
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1999, Liang 2000). The radio halo sizes are typically of (0.5 l)h^o Mpc, and are centered 
on the X-ray emission of the cluster. Their radio halo luminosity at 1.4 GHz, Ji 4, correlates 
strongly with the cluster IC gas temperature, T, and with its X-ray luminosity, Lx, mainly 
produced by thermal bremsstrahlung emission (see Colafrancesco 1999, 2000, Liang et al. 
2000). This result points to a strong correlation between the physical state and origin 
of the relativistic electrons and the thermal content and distribution of the IC gas (see 
Colafrancesco 1999 , 2000 for a more extended discussion). 

While we have a definite evidence for the existence of central radio halos in more than 
20 rich clusters (see Section 6 below), the physical origin of such extended radio halos is 
still a matter of debate. The main difficulty in explaining the radio halo properties arises 
from the combination of their large size (~ 1 Mpc) and the short lifetimes (~ 10^ — 10^ 
years) for synchrotron radiating electrons. Moreover, also the origin and the evolution of 
the IC magnetic field is still unclear (see, e.g., Kronberg 1994 for a review). 

Successful models for the origin of non-thermal, diffuse radio emission from galaxy 
clusters have been proposed in the context of the possible cosmic ray acceleration 
mechanisms in the intra cluster medium (hereafter ICM). Synchrotron emission from a 
population of relativistic electrons (Jaffe 1977, Rephaeli 1977, Roland 1981) accelerated - 
or reaccelerated in situ (Schlikeiser 1987) - as primary cosmic rays (see Sarazin 1999 and 
references therein), or produced in situ as secondary products of pp collisions (Dennison 
1980, Vestrand 1982; see also more recently Blasi and Colafrancesco 1999 and references 
therein) can reproduce the observed spectra of the radio halo emission from Coma (see 
Giovannini and Feretti 1996, Deiss et al. 1997 for an observational review). In all these 
models for the origin of radio halos, one needs a substantial energy supply to maintain the 
brightness of the radio halo for a time comparable with the Hubble time. The possibility 
that such energy budget can be supplied to a radio halo from shocks and turbulence induced 
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by strong cluster mergers (Sarazin 1999, Takizawa 2000) is suggested by the evidence 
that several radio halo clusters (e.g., A1656, A2255, A2163, A2319, A665) are undergoing 
a merger event. The shock waves produced by strong mergers could accelerate primary 
cosmic rays and amplify the IC magnetic field (see Roettiger 1999) to produce substantial 
radio powers. However, there are also some other clusters which still show extended radio 
halos (e.g., A2256, A2319, and other clusters found recently by Giovannini et al. 1999) but 
are either in a pre-merger state or even there is no clear merging event occurring. Moreover, 
there are also a few small-size radio halos found in clusters which have cooling fiows (e.g., 
Perseus, A85, A2218 and maybe A2254). Therefore, one is bound to assume that: i) the 
energy released during the very early phases of the merging process is sufficient to power 
the central radio halo; or ii) the electrons are efficiently reaccelerated by the turbulent gas 
motion originated from IC turbulence and/or the galaxy motion within the cluster (Deiss 
et al. 1997); or Hi) there exists a population of cosmic ray protons accelerated by galaxy 
activity and/or early cluster mergers which are stored for long times (see, Colafrancesco and 
Blasi 1998) in the cluster and produce in situ high-energy secondary electrons providing 
the long-living radio halos (Blasi and Colafrancesco 1999). So, while the energy needed for 
the maintenance of a radio halo can be reasonably accounted for, the origin of the seed 
relativistic electrons is not yet fully clear. 

In this paper we explore the possibility that the diffuse radio halo emission observed 
in many clusters is produced by high-energy electrons resulting from the decay of the 
annihilation products of DM particles, hke WIMP's. Among the WIMP's candidates, 
we consider in particular the lightest neutralino that is predicted in supersymmetric 
extensions of the Standard Model (Haber and Kane 1985, Jungman et al. 1996). The x 
mass is expected to be in the range of tens to hundreds GeV. The detailed interaction 
properties of the neutralino are determined once its physical composition is fixed. Indeed, 
neutralinos are a linear combinations of two Higgsinos (supersymmetric partners of the 
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Higgs bosons) and two gauginos (supersymmetric partners of the neutral gauge bosons). 
Both the neutrahno mass and its physical composition affect in general the rate and the 
final state composition of the XX annihilation process. Recent analysis of accelerator 
constraints, combined with cosmological DM requirements (Ellis et al. 2000), points to: i) 
a lightest neutrahno heavier than about 46 GeV; ii) a lightest neutralino that is either a 
pure gaugino or a heavily mixed gaugino-Higgsino state, since a predominantly Higgsino 
state cannot provide a substantial component of the dark matter. 

According to their physical composition and mass, neutralinos can annihilate either 
into fermion pairs or into vector (or Higgs) boson pairs. In general, for a pure-gaugino x the 
annihilation into fermion pairs is dominant. Due to the chiral structure of the initial state, 
the annihilation cross section inside a halo grows with the square of the final fermion mass 
(see, e.g.. Turner and Wilczek 1990), giving larger rates for heavier fermions. Hence, for 
pure gauginos, a continuum electron spectrum will arise both from the leptonic decays of b 
and c quarks and r leptons, and from the decays of tt's originating from the fragmentation 
of quark pairs. On the other hand, when the x higgsino component is not negligible and 

> Mw, the annihilation into WW (or ZZ) pairs can become dominant. In this case, 
the electron spectrum will arise partly from the leptonic direct W (or Z) decays (giving rise 
to energetic electron of energy ~ M^/2), and partly from the hadronic vector boson decays, 
giving rise to electrons through both the leptonic heavy-quarks decays and tt's decays. 

In this paper, we bypass all the complicancies deriving from a detailed study of the 
supersymmetric parameter dependence of the neutralino annihilation cross sections, by 
adopting the following simplified framework, that will anyway reproduce the crucial features 
of the problem at hand: 

a) we will assume a given value for the neutralino annihilation cross section, which is both 
relevant for the DM closure density and inside the allowed range of the supersymmetric 
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neutralino cross sections for a given and neutralino composition. A different value of 
the cross section can then be straightforwardly rescaled from our results, as we will discuss 
in the paper; 

b) we will study the electron energy spectrum arising from the decays of the neutralino 
annihilation products in the two distinct situations where the dominant annihilation is 
either into fermions or into vector bosons. A realistic case will be, of course, close to either 
these two or to a linear combination of them. 

This simplified scheme will be representative of a wide range of supersymmetric models 
that are relevant for the DM problem. 

The plan of the paper is the following. In Section 2 we discuss the annihilation process 
of neutrahnos in the DM halos of galaxy clusters, and in Section 3 we give an analytical 
approximation for the source spectrum of high-energy electrons produced from the decay 
of the secondary products of XX annihilation. In this Section we also derive the electron 
equilibrium spectrum. In Section 4 we evaluate the radio halo spectrum and its brightness 
profile from synchrotron emission of the relativistic and the level of ICS emission from 
the same population of electrons. In Section 5 we discuss the specific cases of the radio 
halos found in Coma and in the cluster 1E0657-56 and in Section 6 we show how the 
present model is able to reproduce the J1.4 — T relation. We summarize the results of this 
paper and discuss the relevance of the present model for the origin of cluster radio halos in 
Section 7. The relevant physical quantities are calculated using Hq — 50 km s~^ Mpc~^ and 
a flat( Qq — 1) CDM dominated cosmological model. 
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2. Neutralino annihilation in galcixy clusters 

Neutralinos decouple from the primeval plasma when they are no longer relativistic. 
Their present abundance can be calculated by solving the Boltzmann equation for the 
evolution of the density of particle species (see, e.g., Kolb and Turner 1990). 

Assuming a spherical, uniform top-hat model for the collapse of a cluster (see, e.g., 
Peebles 1980), the average DM density within a halo which virializes at redshift z is: 



where the background density is = ^mPc and A(i7oi -2) is the non-linear density contrast 
of the virialized halo (see, e.g., Colafrancesco et al. 1997). Here jg is the gas mass fraction. 
If we redistribute the total mass, M — ^pRp, found within the radius Rp — pre (expressed 
in terms of the cluster core radius Tc), according to a density profile p{r) = pog{r), we find 
that the central total density is given by: 

where I{Rp) = Jq'^ drr'^g{r). In our phenomenological approach, we consider two cases: a 
constant core model and a central cusp profile. 
The constant core model with a density profile: 

2-1 -e/2 



1 + 1- 



(3) 



assumed here to be given from a beta-profile (Cavahere and Pusco-Femiano 1976) with 

^ = 3/3 and [3 = pm.pv'^/kT, gives: 

with p = Rp/vc and I{p,/3) = Jq dxx'^{l + x^)~^^/'^ (here x = r/vc in terms of the cluster 
core radius rc). 
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The central cusp model described by a profile: 



9[r) 



gives: 



Po 



(5) 



(6) 



3 I{P, V, ' 

where I{p,rj,^) — dxx'^~'^ {1 + x)'^~^. With rj — 1 and ^ = 3, the central cusp model 
corresponds to the "universal density profile" which Navarro, Frenk and White (1997; 
hereafter NFW) showed to be a good description of cluster DM halos in N-body simulations 
of CDM hierarchical clustering. Note that, in both models the DM density scales, at large 
radii, as r"^. 

Assuming that the DM density scales hke the total cluster density, a general expression 
for the central DM number density is: 

n3 



np 



(7) 



where the integral / is given in eqs. (4) or (6) for the two DM profiles here considered, and 



n = 4.21 • 10~^ cm"^ Q^/i^ 



r 1 


-1 




LlOO GeV\ 




[ 400 \ 



(8) 



Here we assume that most of the DM consists of neutralinos of mass M^. 
The annihilation rate of neutralinos in a DM halo is 

R^n^{r){aV)A, 



(9) 



where %(r) = n^fig{r) is the neutralino number density and {aV)A is the XX annihilation 
cross section averaged over a thermal velocity distribution at freeze-out temperature. 
Although, as anticipated, the xX annihilation cross section is a nontrivial function of the 
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mass and physical composition of the neutrahno, to our purpose it suffices to recall that the 
Xrelic density is approximately given by ( Jungman et al. 1996) : 

^ , 9 3 X 10-27 cm3 s-i , , 

"^"^ — m-. — ■ 

Hence, for the XX annihilation, we will assume a total x cross section of 

{aV)A ^ 10"^^ cm^ s"^ (11) 

to be consistent with the value flm ~ 0.3 derived from clusters of galaxies and large scale 
structure constraints (see, e.g., Bahcall 1999). Detailed studies of the relic neutrahno 
annihilation (Edsjo 1997) show that the above value is well inside the allowed range 
predicted in supersymmetric theories for a wide choice of masses and physical compositions 
of neutrahnos that can be relevant as cold DM candidates. Enhancing the x annihilation 
rate will have on our results the simple effect of rescaling the final electron spectra by the 
same enhancement factor 

For values of h in the range 0.5 ^ /i ^ 1, we can infer that for a fiat, vacuum dominated 
CDM universe with Q\ in the range 0.5 -r- 0.7 and ~ 0.3 -r- 0.5, the quantity Q^h"^ takes 
values in the range ~ 0.075 -j- 0.5, which would fix the annihilation cross section to within a 
factor less than ten. 



3. Production and equilibrium electron spectra 

3.1. The source spectrum 

Neutralinos which annihilate in a DM halo produce quarks, leptons, vector bosons and 
Higgs bosons, depending on their mass and physical composition (see, e.g., Edsjo 1997). 
Following the discussion in Section 1, monochromatic electrons (with energy about M^), 
coming from the direct channel xX ~^ ^e, are in general much suppressed (Turner and 
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Wilczek 1990). Electrons are then produced from the decay of the final heavy fermions and 
bosons. 

The different composition of the XX annihilation final state will in general affect the 
form of the final electron spectrum. To our purpose, it will be sufficient to consider two 
somewhat extreme cases: 1) a pure-gaugino annihilation, which yields mainly fermion pair 
direct production xX ~^ f f i with cross sections scaling as M| for different fermions (Turner 
and Wilczek 1990); 2) a mixed gaugino-higgsino state in the case of a dominant annihilation 
into W (and Z) vector bosons, xX ~^ WW{ZZ). A real situation will be mostly reproduced 
by either of the above two cases, or by a linear combination of the two. 

The positron spectrum arising from the xX annihilation has been derived by various 
authors (Silk and Srednicki 1984, Rudaz and Stecker 1988, Ellis et al. 1989, Stecker and 
Tylka 1989, Turner and Wilczek 1990, Kamionkowski and Turner 1991, Baltz and Edsjo 
1998). Here, we will adopt the approaches by Rudaz and Stecker (1988; hereafter RS) and 
by Kamionkowski and Turner (1991; hereafter KT), that gave the analytical approximations 
of the positron source functions for models in which neutralinos annihilate mainly into 
fermions and vector bosons, respectively. Here we report the results of our calculations 
while the details of these two approaches are presented in the Appendix. 

The total source spectra in the two cases considered are shown in Fig.l for 
= lOOGeV. The heavy solid curve is the total spectrum for the model in which fermions 
dominate the annihilation. We also show the different contributions to the total source 
spectra from fermions, i.e. the source spectrum for first-generation prompt electrons (PI), 
second-generation prompt electrons (P2) and secondary electrons produced in the decay of 
charged pions (tt). The total source spectrum is rather smooth and, when approximated 
by a single power-law, it has an overall slope Qe ~ E~^-^ in the interesting energy range 
0.02 < E/M^ < 0.7. 
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The source spectra for the model in which gauge bosons dominate the annihilation 
is given by the light solid curve. For this case we show the contributions from the decay 
t"^ —>■ e^; — > c — > (labelled as KTi) and from the decays of charged pions 
produced in hadronic decays of the (labelled with KTtt). Note that we make this 
analysis restricting only to the 1^1^ channel, since the ZZ contribution gives rise to a 
qualitatively similar spectrum. 

Two prominent spectral features (bumps) around E ~ My./2 (arising from the direct 
decay W ev) and E ~ are shown by the model in which gauge bosons dominate 

the annihilation, in comparison with the smoother spectral shape of the fermion dominated 
annihilation. These spectral features are more prominent in Fig.2 where we plot the 
quantity E^Qe{E,r). From a closer inspection of Figs.l and 2, one can see that the 
main difference between the two models stands in the height and width of the bumps at 
energies below M^/10, where the pion-produced electrons dominate the source spectrum, 
and in the high-energy tail at £■ ^ O-^-Z^c where the PI contribution dominates. In both 
models, electrons produced by tertiary decays of gauge bosons arc neglected because their 
contribution is subdominant with respect to the pion decay distribution. 

In the calculation of the radio halo emission (sec Section 4 below) the source spectra of 
eq.(Al) and A(18) will be multiplied by a factor 2 to take into account the contribution of 
both electrons and positrons. 



3.2. The electron equilibrium spectrum 



The time evolution of the electron spectrum is given by the transport equation: 

dne{E,r) d 



dt 



dE 



ne(E,r)b{E) 



^Qe(E,r) 



(12) 
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where ne{E,r) is the equilibrium spectrum at distance r from the cluster center for the 
electrons with energy E. The source electron spectrum rapidly reaches its equilibrium 
configuration mainly due to synchrotron and Inverse Compton Scattering (hereafter ICS) 
losses at energies E ^ 150 MeV and to Coulomb losses at smaller energies (Blasi and 
Colafrancesco 1999). Since these energy losses are efficient in the ICM and DM annihilation 
continuously refills the electron spectrum, the population of high-energy electrons can be 
described by a stationary transport equation {dne/dt ~ 0) 



_d_ 
dE 



neiE,r)b{E) 



Qe{E,r) (13) 



from which the equilibrium spectrum can be calculated. Here, the function b{E) gives the 
energy loss per unit time at energy E 

where bo{B^) = (2.5 • 10-^^ + 2.54-10-^^52) and bcoui = 7 x lO-^S(r) (if be is given in units 
of GeV/s). In the expression for bcoui, the IC gas density, n(r), is given in units of cm~^. 

The equilibrium spectrum calculated combining the source spectrum Qe{E, r) and 
the energy losses b{E) is shown in Fig. 3 for the two cases here considered. Equilibrium 
spectra are evaluated separately for the different components of the source spectra and 
then summed together. This spectrum is evaluated within the cluster core for a constant 
core density (see eq.2) with a reference value of the central density n^_o — 1 cm~^ and 
{(jV)a = 10"^^ cm^ s"^. This equilibrium spectrum can be easily calculated for other cluster 
configurations by scaling the parameters of the cluster density profile. The equilibrium 
spectrum for the case in which neutralino annihilation is dominated by fermions, is smooth 
enough that it can be reasonably approximated by a power-law, ne{E) ~ E~p, down to 
E ~ M^/30 with a slope p ~ 2.9. We will make use of this convenient approximation in the 
discussion of our results. 
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3.3. Comparison between approximated and Monte Carlo results 

The source spectra in Fig. 1 arise from an analytic approximation (see the Appendix) 
of the exact shape of the electron spectrum that tries to cope with the details of the quarks 
and leptons decays and of the hadronization of the decay products. Detailed spectra can 
also be obtained by using state-of-the-art Monte Carlo simulations, although the analytical 
approximations used here are able to recover the relevant aspects of more detailed studies. 

There are, however, some differences among the equilibrium spectra obtained from 
different studies using Monte Carlo techniques. We can compare here our analytical 
approximations with the results of the Monte Carlo simulations run by Baltz and Edsjo 
(1998), and Golubkov and Koplinoch (1998). Note that the latter study considers massive 
neutrinos annihilation, that, as far as the electron spectra shapes are concerned, should be 
equivalent to the neutralino annihilation at fixed = rriv, and fixed ff /WW annihilation 
dominance. One can note that for the case we consider here, i.e. = 100 GeV, the total 
equilibrium electron spectrum from the simulations of Golubkov and Konoplich (1998) 
agrees quite well with the spectrum shown in Fig. 3 for the case in which neutralinos 
annihilate mainly into fermions. On the other hand, the case in which the annihilation 
is dominated by vector bosons (KT) is quite well reproduced by the simulations in Baltz 
and Edsjo (1998), although there a slightly different value of (i.e., ~ 130 GeV) is 
assumed. 

However, beyond the limits of the present discussion, it is worth noticing that there 
are still differences in the electron spectra so far pubhshed by different authors using Monte 
Carlo techniques. In view of this, we conservatively assume in our calculations a reasonable 
uncertainty of a factor ~ 2 in the amplitude of the electron spectrum, and of ~ 10% in 
its slope. We will consider more detailed predictions derived from available Monte Carlo 
simulations in a forthcoming paper. 
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4. The diffuse radio halo emission 

The calculation of the radio emissivity per unit volume is performed here in the 
simplified assumption that electrons with energy E radiate at a fixed frequency given by : 

i/Ri3.7MHz5^(^^y . (15) 

This approximation introduces negligible errors in the final result and has the advantage of 
making it of immediate physical interpretation. 

The radio emissivity at frequency v and at distance r from the cluster center can be 
calculated as 

jM=ne{E,r)[ — \ —, (16) 

\ / syn 

where {dE/dt)syn is given in eq.(14). The radio halo luminosity is obtained by integration 
of j(i/, r) over the cluster volume yielding: 

[Rhalo 

J(u) =A7r dr r^j(iy,r) , (17) 

Jo 

where Rhaio is the size of the cluster radio halo. The observed radio halo fiux is then 

J{u) 



where Dl{z,0,o) is the cluster luminosity distance (see, e.g., Weinberg 1972). For the sake 
of illustration, we give the explicit form of J{i') in the case of a power-law source spectrum, 

Q{E,r) = Qq{E / Eq)^% where Qo = n'^{r){aV)A{kM^)^^ and Eq is a normalization energy, 
say 1 GeV. Using eqs. (13), (16) and (17), we can write 

W^) y'^o) ^ ^ 

where = s/2, E, = 16.44 GeV, J = /^f drr^g^{r), ho{B^) = (2.5 ■ IQ-^^ + 2.54 ■ 10-^^3^ 
and bo^syniBfx) = 2.54 • 10~^^B^. For ^ 3.14, the radio luminosity scales in eq.(18) as 
J(u) (X nl^o{(^V)A{kM^)-^B^''+'^u-''''. 
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In the calculation of the radio halo luminosity we consider both the case of a uniform 
IC magnetic field and a case in which the magnetic field has a radial profile which declines 
with the distance from the cluster center, as indicated by recent numerical simulations (see, 
e.g., Priaca and Goncalves 1998). Specifically, we considered a magnetic field density profile 



B = BM 



^ ^ 2-| -w/2 



with values i/; ~ 0.5 -j- 1. Here, for simplicity, we assume that the spatial structure of the 
magnetic field is given by a constant core profile with lenght scale Tcb- 



4.1. The radio halo spectrum and spatial distribution 

The total radio halo flux predicted from our model has generally a steep spectrum 
with average slope ~ 1.2 — 1.8, depending on the neutralino mass and composition 
and on the shape of the DM density and magnetic field profiles. We will discuss more 
specifically this point in Section 5 below. Here we want to emphasize that for a fixed value 
of the magnetic field the neutralino mass fixes the maximum frequency at which we can 
observe synchrotron emission from the radio halo, 

i/„„,~3.7MHzB^(^^j (20) 

For = 100 GeV and k ^ 0.7 (see Appendix), the maximum frequency at which the 
radio halo spectrum is observable is I'max ~ 18.1 GHz x5^. From this fact, we can use the 
available data on the radio halos of Coma and 1E0657-56 (see Figures 3 and 4) to set a 
lower limit to the neutralino mass. Specifically, the maximum frequency, i'max,obs, at which 
the radio halo spectrum is observed sets a lower limit 

M, > 1^ GeV (^^) '^'s;V2 (21) 

which gives > 74.3 GeV B~^/'^, for Umax,obs — 10 GHz. Furthermore, the present 
theory predicts a sharp cutoff in the radio halo spectrum at Vmax given by eq.(20); so. 
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high frequency observations can both test the neutrahno model and measure the quantity 
B^{kM^)'^ from the possible detection of the high-i/ sharp cutoff in the radio halo spectrum. 

The brightness profile of the radio halo is obtained by integration of the emissivity, 
_7(z/, r), along the line of sight. It scales, for a power-law radio halo spectrum, as 

Sradioo^ J den\r)-B'+''^{r) (22) 

and declines at large radii as 7--K+«'{i+"r)]+i/2 ^j^g radio halo brightness profile resembles, 
for both the DM profiles here considered, the X-ray brightness profile of the cluster which 
scales as 

Sx-rayo^ J M n\r) ■ T^''' {t) . (23) 

In particular, for a uniform magnetic field, B — const, and for a isothermal cluster, T — 
const., the brightness S radio behaves exactly as the X-ray brightness S x-ray ^ radii. 
Small deviations from this behaviour can be attributed to radial variations in B{r) and/or 
in T{r). 

The radio halo brightness profile depends on the assumed DM profile (see eqs.19 and 
22). Specifically, at small radii the constant core model still gives a radio brightness profile 
that resembles the X-ray emission profile while the central cusp model tends to give a radio 
brightness which is more peaked toward the center. Thus, the model predicts a strong radio 
emissivity at the cluster center under the assumption of a central cusp (or NFW) profile. 
However, neutralino annihilation can also cause a softening of the central DM cusp due to 
the high interaction rate of neutralinos in the central region of the clusters with respect to 
the outer parts (see Kaplinghat, Knox and Turner 2000). This countereffect weakens the 
aforementioned problem. 
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4.2. Inverse Compton Scattering emission 



The relativistic electrons which are responsible for the radio halo emission also 
emit X-rays and UV photons through Inverse Compton Scattering (ICS) off the Cosmic 
Microwave Background photons. As in the case of the synchrotron emission, also for ICS 
we can adopt the approximation that electrons radiate at a single energy, given by 



Electrons with energy in excess of a few GeV radiate in the hard X-ray range, while 
electrons with energy smaller than ~ 400 MeV produce soft X-rays and UV photons. The 
non-thermal X-ray/UV emissivity at distance r from the cluster center is evaluated as 



In complete analogy with the case of the radio emission, the integrated non-thermal X-ray 
luminosity, $x(-E'x) is 



The predicted X-ray fluxes from ICS emission due to the electrons produced in XX 
annihilation are lower than those required to explain the hard X-ray excesses currently 
observed in Coma, A2256 and A2199, and do not contribute substantially to the EUV and 
hard X-ray emission excesses of these clusters. 

5. Application to observed radio-halo clusters: Coma and 1E0657-56 

Our theory for the origin of radio halos in clusters can reproduce successfully the 
spectra of the radio halo emission observed from Coma (Giovannini et al. 1993, Deiss et al. 
1997) and 1E0657-56 (Liang et al. 2000). Moreover, these data in turn can set interesting 
constraints to the mass and composition of the neutralino. 



Ex ^2.7 keV E^GeV). 



(24) 




(25) 




(26) 
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For any cluster, we choose a DM central density, n^fl, which is given by eq. (7) 
evaluated with the parameters appropriate to each cluster and, in the case of 1E0657-56, 
with parameters approriate to each reagion in which the radio halo spectrum has been 
observed. 

In the case of Coma we choose the following parameters: Tc = 0-4:h^Q Mpc, 
P — 0.76, riy^fi — 3 • 10"^ cm~^. The radio halo spectrum has been integrated out to a 
radius Rhaio — ^-^h^o Mpc (see, e.g., Giovannini et al. 1993). In the model in which 
neutralino annihilation is dominated by fermions {xx ~^ f f )i ^ uniform magnetic field 
^uniform ~ 1-3 /iG is needed to fit the spectrum under the assumption of a constant core 
profile. In the case of a dechning magnetic field B — So[l + (r/rc,B)^]~°'^, a value Bq — S /iG 
fits better the data, under the same model assumptions. The resulting spectra are shown 
in Fig. 4. In this figure we also show the spectra derived from the power-law approximation 
of the equilibrium spectrum shown in Fig.3. A radially decreasing magnetic field yields a 
slightly steeper spectral slope, as it should be expected from eq.(17). 

Under the same assumptions of a constant core profile and radially decreasing 
magnetic field, the model in which neutralino annihilation is dominated by gauge bosons 
{xx ~^ WW) does not fit the Coma radio halo (see Fig.4). Within the limits of the 
analytical approximations used in this paper, the Coma data exclude this model at more 
than 3 standard deviations (see Fig.5). In this figure we also show the range of uncertainties 
in the final radio halo spectra calculated considering an uncertainty of a factor ±2 in the 
overall normalization of the source spectra and ±10% uncertainty on their slopes, for the 
two extreme annihilation model here considered. 

In the case of the distant cluster 1E0657-56 at z = 0.269, we choose different parameters 
according to the two different regions in which the spectrum has beem measured (we use 
the core radii and /3 parameters given in Liang et al. 2000). We assume that neutrahno 
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annihilation is dominated by fermions (xx — > //) and consider a constant core DM density 
profile. For the outer region we choose Tc = O.SSh'^Q Mpc, P = 0.7, n^^ = 9 • 10~^ cm~^. 
The radio halo spectrum has been integrated out to a radius Rhaio = '^h^o Mpc (Liang et al. 
2000). A uniform magnetic field Buniform ~ 2. /iG is needed to fit the spectrum of the outer 
region. In the case of a declining magnetic field B — Bo[l + (r/rc,B)^]~°"^, a value Bq — 100 
IJ,G can fit the data. For the inner region we choose = O.OSh'^Q Mpc, (3 = 0.49 and the 
same central density n^fi — 9 • 10~^ cm~^. The radio halo spectrum has been integrated out 
to a radius Rhaio — O.B/igo^ Mpc. A uniform magnetic field B Ki 9 /iG is needed to fit the 
spectrum of this region. In the case of a declining magnetic field B = Bq[1 + (r/rc,B)^]~°'^, 
a value Bq = 100 /iG fits the data. The results are shown in Fig. 7. As in the case of Coma, 
the extreme case of neutralino annihilation dominated by gauge bosons (xx ~^ WW) is 
inconsistent with the radio halo spectrum of this cluster. 

Note that the radio-halo spectra shown in Figures 4-7 are calculated assuming a 
constant core density profile for the DM component. A central cusp model gives a higher 
central density and increases the x annihilation rate in the central region of the cluster 
which in turn produces more high-energy electrons per unit volume. This increase in the 
electron density must be compensated by a decrease in the magnetic field value needed to 
fit the radio halo spectrum. In fact, assuming a NFW density profile for Coma, a uniform 
magnetic field Buniform = 0.45/xG is required to fit the data. Alternatively, a central value 
of the magnetic field, Bq — 1.4/xG is able to fit the spectrum assuming a B radial profile 
as given by eq. (19). We show the difference between the spectra obtained either with a 
beta-profile or with a NFW profile in Figure 8. 
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6. The Ji.4 — T relation for radio-halo clusters 

Another interesting property of radio halo clusters is the steep correlation existing 
between the monochromatic radio halo luminosity observed at 1.4 GHz, J1.4, and the IC gas 
temperature: J1.4 ~ with q = 6.25 (range g = 4.1 -^ 12.5 at 90% confidence level) found 
by Colafrancesco (1999, 2000) and confirmed by Liang et al. (2000). Radio luminosity 
data are taken from Feretti (1999), Giovannini et al. (1999), Liang et al. (1999) and Owen 
et al. (1999). Here, we use the temperature data from Arnaud and Evrard (1999), bur a 
similar result obtains using temperature data from Wu, Xue and Fang (1999) or David et 
al. (1993). Such a steep correlation (see Fig. 9) can be reproduced in models of radio halos 
powered by acceleration of cosmic rays triggered either by strong merging events or by 
enhanced galaxy activity, but requires specific conditions for the structure of the magnetic 
field (see Colafrancesco 1999, 2000 for a discussion). 

In the model we present here the radio halo luminosity (see eqs. 17-18) scales as 

J{u) oc nl^orl{aV)AB^+'''-iy-'''- (27) 

where ar = dlnJ{u) / dlnu is the effective slope of the radio halo spectrum. Thus, 
the monochromatic radio halo luminosity correlates naturally with the cluster X-ray 
bremsstrahlung luminosity, Lx oc n^fgr^T^/'^ (here no is the central IC gas density and 
fg = Mgas/M is the cluster baryon fraction), to give 

J1.4 cx Lxf;^T-^'\aV)AB^^'''^ . (28) 

Using the observed correlation, Lx ~ with 6 ~ 3 (see, e.g., Arnaud and Evrard 1999, 
Wu et al. 1999, David et al. 1993) and assuming that fg = const, as indicated by the 
available cluster data (Jones and Forman 1999, Mohr 2000), we derive 

JiA ~ r(^''-^)/2 ^^v)aB^+"^ (29) 
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Under the condition - assumed in this paper - of hydrostatic equihbrium of the cluster 
material with the overall potential wells, the relation V ~ T^/^ holds, and we finally obtain: 

Ji.4 ~ . (30) 

If every cluster has the same value of the magnetic field (universality condition), a 
correlation J1.4 ~ obtains. This is too flat compared to the observed correlation. 
However, the universality condition B — const is rather unlikely. In fact, we expect that B 
increases with increasing T either under the condition of energy equipartition between the 
IC gas and the magnetic field (which yields B ~ T^/'^) or even under the condition of a 5 
field which is frozen- in the ICM (i.e., B n' with 7 ~ 2/3, which yields B - T^''-^)^). 

In conclusion, it is reasonable to expect a quite steep correlation, Jia ~ 7'4.i5±o.25 
or J1.4 ~ 2^4.4o±o.26^ equipartition or frozen- in magnetic field. Such a correlation is 

consistent (within the uncertainties) with the best fit result given above and can reproduce 
the overall distribution of clusters in the J1.4— T plane. Moreover, since T oc {l+z)/S^/^{ytQ, z) 
in CDM models of structure formation (see, e.g., Colafrancesco et al. 1997), the radio halo 
power is expected to evolve strongly with redshift J1.4 ~ [(1 -|-2;)A^/^(Oo, Thus, already 
at 2; ~ 0.2 its normalization increases by a factor ~ 2 (we use here = 4 as a reference 
value). The Jia — T correlation predicted in the present model at 2; = and at 2; = 0.25 is 
shown in Fig. 9. The large steepness q ~ 6.25 deduced by a formal power-law fit to the data 
can be due to a superposition of the different correlations of radio halo clusters observed in 
the range 2; 0.3. 

7. Discussion and conclusions 

In this paper we have shown that the basic properties (i.e., the spectrum, the 
surface brightness profile and the J1.4 — T correlation) of radio halos in galaxy clusters 
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can be reasonably fitted by a model in which the high-energy electrons, giving rise to 
the synchrotron emission, arise from the decay of secondary products of the neutralino 
annihilation in the DM halos of galaxy clusters. The slope and the frequency extension 
of the spectrum can set interesting constraints on the neutralino mass and composition. 
In fact, the highest observed frequency of the radio halo spectrum set a lower hmit to the 
neutralino mass: 



In particular, the highest frequency {u ~ 4.85 GHz) at which the Coma spectrum has been 
observed requires > 54.6 GeV B~^/'^. A more stringent limit, > 70.5 GeV B~^/'^, 
is obtained from the cluster 1E0657-56. Radio halo observations at frequencies larger 
than 5 -i- 10 GHz are not yet available, but are important to test the present model. In 
fact, the detection of a high frequency cutoff, Ucut, in the radio halo spectra gives an 
upper limit on the neutralino mass. The high-i/ cutoff predicted in our model, has the 
same value, Ucut — 3.7 MHzS^[A;(M^/Gey)]^, for all radio halos and depends only on the 
neutralino mass and the IC magnetic field. For a fixed neutralino mass and composition, 
the observations of this high-z/ cutoff can yield a direct estimate of the cluster magnetic 



In the same framework, the slope of the radio halo spectrum gives an indication on the 
neutralino physical composition. The radio halo spectrum of the two clusters considered 
here are well fitted by a model in which the neutralinos behave like pure gauginos, and 
annihilate mainly into fermions. On the other hand, when the neutralino annihilation is 
dominated by vector bosons (implying a non-negligible higgsino component), the electron 
source spectrum is too steep and shows two unobserved bumps at low and high energies 
(see Figs.l and 2). This feature is also evident in the final radio halo spectrum (see Fig. 5). 
Indeed, Figs. 5 and 6 show how far the neutralino model dominated by vector bosons is 
from reproducing the data of the radio halo spectrum of Coma and 1E0657-56. 




field. 
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After setting the neutralino mass and composition, the slope of the radio halo spectrum 
still depends on the DM density profile and on the magnetic-field radial profile, that 
influence the annihilation rate and the synchrotron emission power, respectively. The 
available data on Coma indicate (in the case of a constant core density profile) that a 
radially decreasing magnetic field is favoured to fit the correct steepness of the spectrum 
at high frequencies. This constraint is tighter for Coma, where an upper limit of 52 mJy 
at u — 4.85 GHz has been measured (see, e.g., Deiss et al. 1997), than for 1E0657-56, 
where a rather constant slope is observed up to ~ 10 GHz. For a constant core density 
profile, the spectra of Coma and 1E0657-56 are consistent either with uniform IC magnetic 
fields Buniform ~ 1.3 -j- 2 /iG or central values Bq 8 -r 100 iiG for a magnetic field profile 
B{r) ~ at large distances from the cluster center. Using a central cusp (NFW) density 
profile increases slightly the slopes of the radio halo spectrum, reducing consequently the 
need for a magnetic field strongly decreasing with the distance from the cluster center. The 
central cusp profile requires also IC magnetic field amplitudes lower by a factor ~ 3 -^ 5 to 
fit the radio halo spectrum. 

The present model also predicts a mild steepening of the radio halo spectrum with 
increasing distance from the cluster center. This effect is mainly due both to the lower 
DM densities (which decrease the neutrahno annihilation rate) and to a decrease of the 
IC magnetic field (which reduces the emitted powers via synchrotron emission) at large 
distances from the cluster center. The model presented here can fit reasonably well the 
radio halo spectra of the cluster 1E0657-56 in the inner and outer regions at which it has 
been measured by Liang et al. (2000). It also gives values for the radio spectral index of 
Coma ar ~ 0.9 -j- 1.6, from the cluster center out to ~ 1 Mpc. These values are consistent 
(within the given uncertainties) with the increase of the radio spectral index of Coma given 
by Giovannini et al. (1993). 
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It is also noticeable that the spatial extension of the radio halo emission which is 
predicted in the present model is very similar to that of the X-ray emission from the IC gas, 
a property which has been emphasized in several observational works (see, e.g., Feretti 1999, 
Liang et al. 2000). Changes of radio halo/X-ray surface brightness ratio, Sradio/ S x-ray, 
with radius can be accounted for by radial variations of the IC magnetic field and/or the 
IC gas temperature. 

The analysis presented in this paper is made for a value of the neutralino annihilation 
cross-sections {(tV)a ~ 10^^^ cm^'s^^. This value is well inside the allowed range for a 
neutralino relevant as a DM candidate, and can be rescaled in a straightforward way from 
our results, for different assumptions on {aV)A- Because the radio halo spectrum power 
(see eqs. 17-18) scales hnearly with (aV) An'^^B^^"'' (here = dlnF{y) / dlnv is the effective 
spectral slope) , higher (lower) values of the annihilation cross section imply lower (higher) 
DM densities, and require a lower (higher) IC magnetic field. 

Another interesting aspect of the present model for the origin of cluster radio halos 
is that it can reproduce fairly well the observed correlation between the monochromatic 
radio halo power, J1.4, and the IC gas temperature T. The steep J1.4 — T relation shown 
by the available radio-halo clusters can be reproduced in our model as a superposition of 
evolutionary effects of the correlation J1.4 ~ T^, with q ~ 4.2 -^ 4.4, which results from 
the dependence of the electron spectrum, ne{E,r) oc n^{r){aV) a, from the DM density, 
n^(r), and annihilation cross-section, {(jV)a, with the further assumption of an energy 
equipartition between the IC gas and the IC magnetic field. 

The ICS emission from the same population of relativistic electrons produces 
unavoidably also fluxes of UV and X-ray emission that, however, are not so intense to 
reproduce the emission excesses in the EUV (Lieu et al. 1996) and in the hard X-rays 
(Fusco-Femiano et al. 1999, Rephaeli, Gruber and Blanco 1999) observed in Coma. 
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This fact should be not considered as a problem for the model we have worked out here 
since other alternative explanations have been proposed to fit the EUV and hard X-ray 
emission excesses in the framework of thermal (see Antonuccio-Delogu et al. 2000) and/or 
suprathermal phenomena (Ensshn et al. 1999-2000, Dogiel 1999, Sarazin and Kempner 
1999, Blasi et al. 2000), respectively . 

Another remarkable feature of the present model, is that neutralino annihilation can 
also give rise to gamma rays with continuum fluxes which are overwhelmingly due to 
7r° — > 7 -|- 7 decays. The continuum gamma-ray spectrum is given by: 

f\{E,) = 2 / dE^{El - ml)-'/%f{E^) (31) 

(see,e.g., Rudaz and Stecker 1988) where the quantity C7r/(-E'7r) is derived in equation (A12) 
of the Appendix. We recall here that E^^ is the pion energy, Ei = E^ + ml/iE^ and E-y is 
the energy of the gamma-ray photon. The resulting gamma-ray luminosity is 

L, = A7rnl,{aV)A j"" drr'g\r)f,{E,) (32) 

J 

and yields, for Coma, a gamma-ray flux F^{> lOOMeF) = 3 • 10~^ cm'''^s~^ (for 

= 100 GeV; see Fig. 10), which stays below the EGRET upper limit for Coma, 
pComa(^^ IQQMeV) = 4 • 10"^ cm-^s"! (Sreekumar et al. 1996). 

Note that a correlation oc ^■'^^^ (here b is the exponent of the X-ray luminosity 
- temperature correlation, Lx ~ T^) is expected in this model with fg — const, at variance 
with the correlation L-y ~ expected in models of cosmic ray (pp) interaction in the 
ICM, (Colafrancesco and Blasi 1998). The different slope of the — Lx relation will 
provide a way to disentangle between these two mechanisms for the production of gamma 
rays in galaxy clusters. 

The gamma ray emission produced by XX interaction extends up to energies 
corresponding to the neutralino mass (see Fig. 11). One way to disentangle a particular 
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model for the cluster gamma-ray emission is provided by the future high-energy gamma-ray 
experiments. They will allow to observe galaxy clusters with a photon statistics sufficient to 
disentangle the gamma-ray spectra produced by the 7r° ^ 7 -|- 7 electromagnetic decay, as 
predicted either in DM annihilation models {xX ~^ ff: WW — > 7r°) or in secondary electron 
models {pp — > n^; see, e.g., Colafrancesco and Blasi 1998), from the gamma-ray emission 
produced by the bremsstrahlung of primary cosmic ray electrons. These observational 
capabilities will be available, however, in the next future. We will discuss more specifically 
this issue in a forthcoming paper. 

Finally, we want to discuss a number of additional stringent predictions through which 
the model discussed here can be tested. 

(i) According to the fact that DM is present in all large scale structures, we should observe 
radio halos in every cluster, which is not actually observed. In the framework of this paper, 
we have shown that radio halos can be fitted with sensitively high magnetic fields (with 
central values ~ 5 100 //G). Then, only clusters which have such high magnetic fields 
can show a bright radio halo. All the other clusters are expected to have faint radio halos 
that can brighten up when there is some effect that raises the magnetic field amplitude, for 
example, a merging event (see the numerical simulations of K. Roettiger. 1999). 

(ii) A possible problem could be given by the fact that cooling flow clusters, which have 
usually high magnetic field in their central cooling regions, do not show strong evidence 
for extended radio halos. However, if one assumes that the magnetic field in cooling fiow 
clusters is quite peaked near the cluster center and decreases rapidly towards the outskirts, 
the radio halo emission predicted in the present model is a factor ~ 10^ -7- 10^ lower than 
that of a non-cooling fiow cluster with the same mass [these estimates are obtained for 

a Coma-like cluster using the magnetic field profile given in eq.(19) with parameters, 
fc,B = O.l/i^o^ Mpc and w = 0.5 2]. Nonetheless, it remains true that our model does 
predict that cooling fiow cluster should possess small-size, low-luminosity and low-surface 



-29- 



brightness radio halos. 

(in) Under the assumption that all clusters show a universal DM profile, the radio halo 
spectrum is basically the same unless the magnetic field configuration of the cluster is very 
peculiar. For example, steep radio halo spectra can be obtained mainly due to the presence 
of magnetic fields that decrease strongly from the cluster center. 

Little is known about the presence and structure of the IC magnetic fields. Faraday rotation 
measurements (Vallee et al. 1986, 1987) give a lower estimate of the magnetic field in small 
scale regions of the cluster (Dolag et al. 1999), and we could have to deal with large scale IC 
magnetic fields that are lower than those {Bo ~ 5 -j- 100 fiG, Buniform ~ 1 ^ 3 jiG) required 
by the present model to explain radio halos. However, even though the magnetic field are 
not so strong in clusters, the present model certainly yields a natural explanation for the 
origin of the seed high-energy electrons, which are a necessary input for any model of radio 
halo (and relic) formation. The diffuse radio emission due to such a population of seed 
electrons could then be boosted by the amplification of the IC magnetic field subsequent to 
a strong cluster merger (Roettiger 1999), or the seed electrons could be reaccelerated by 
intracluster turbulence (Deiss et al. 1997, Eilek and Weatherall 1999, Brunetti et al. 1999). 

While at the moment, we observe radio halos in more than 20 clusters at different 
redshifts, the definite test for the theory of the radio halo origin proposed here is committed 
to obtain a larger, unbiased survey of galaxy clusters through high radio sensitivity 
observations. This search should be complemented with the search for gamma-ray emission 
from galaxy clusters whose predicted intensity is matched by the sensitivities of the 
next generation gamma-ray experiments (GLAST, AGILE, MAGIC, VERITAS, ARGO, 
STACEE). 

In conclusion, we want to emphasize that the astrophysical expectations from the XX 
annihilation are consistent, at the moment, with the constraints set by all the available 
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observations on clusters containing radio halos. The astrophysical and fundamental physics 
requirements on the model discussed here are stringent, but still well allowed. It is also 
appealing, in these respects, to expect that some astrophysical features of galaxy clusters 
might give information on the fundamental properties of the DM particles. 
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A. APPENDIX. The production spectrum of high-energy electrons 

In this Appendix we will derive an analytical approximation for the source spectrum of 
electrons resulting from the decay products of XX annihilation. Here we will refer mainly to 
the analytical approaches of Rudaz and Stecker (1988; hereafter RS) and of Kamionkowski 
and Turner (1991; hereafter KT) who gave analytical approximations of the positron source 
functions for models in which neutralinos annihilate mainly into fermions {xX ~^ ff) o^" 
gauge bosons {xX ~^ WW) , respectively. In the following, we derive the electron source 
spectra for generic values of the neutralino number density, n-^{r), and annihilation cross 
section, {aV)A- 



A.l. Fermion dominated annihilation 

Following RS, we consider three main sources of from xX annihilation: (PI) first 
generation, prompt electrons with a continuum spectrum; (P2) second generation, prompt 
electrons; (tt) electrons produced from the decay of tt^. The electron source spectrum, 
Q{E, r), can be written as the sum of these three components 

Qff(E,r)^Qpi(E,r) + Qp2(E,r) + Q^(E,r) . (Al) 

For a neutralino with mass = 100 GeV, the first generation source spectrum takes 
contributions mainly from 6 — > e, c ^ e, ^ (see RS). Taking into account the different 

energy distributions of the previous decays, the first generation electron spectrum is found 
to be well approximated by 

Qpi{E,r)^nl{aV)ACifi{E) (A2) 
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with Ci ^ 0.6/{kM^) and k ^ 0.7. 
The function fi{E) is given by 



where 



Mx) = ^fr{x) + 0.13 X + 0.1 X '^Mx) (A3) 



Mx) = Mx) = l + lx'- (A4) 



and 

fdx) = 2(^l + 2x^ - 3x^^ (A5) 

in terms of the adimensional quantity x — E/M^. In eq.(A3) we replace the values assumed 
by RS for the leptonic decay branching ratios with the present values taken from the 
Particle Data Group (Groom et al. 2000). Neglecting the differences between fr, fb and fc, 
the source function Qpi - as noticed by RS - can be approximated by 

Qpi(E, r) ^ 0.nnl{aV)A ■ f(E) (A6) 

where 

m = j^JikM^ - E) (A7) 

where the function 9{kM^ — E) is the Heaviside function. The source function Qpi{E,r) is 
basically constant up to E ^ kM^ and then drops to zero for higher energies (see Fig.l). 



The second generation source spectrum (arising from the leptonic decays in the final states 
of the first generation decays) is found to be well approximated by: 



Qp2{E,r)=nl{aV)AC2f2{E) 



(A8) 



with 



f2iE) 



km. 



-In 



. E 



(A9) 
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where, again, k ~ 0.7 and 

17 3 27 

C2 ~ + 0.13 X — + 0.1 X 2 X — . (AlO) 

This source function goes rapidly to zero for E — > fc^M^ ^ -^x/^ (see Fig.l). Again, in 
eq.(AlO) we replace the values assumed by RS for the leptonic decay branching ratios with 
the present values taken from the Particle Data Group (Groom et al. 2000). 



Following RS, we find that the source spectrum from tt^ decay is found to be well 
approximated by: 

Q^{E,T)=nl{aV)AU.{E) (All) 

with 

U{E) = e-°-6«^ + 0.115e-°-2^«^ (A12) 

and Ctt - 1.74(100 GeV/M^). 

The three source functions for PI, P2 and tt contributions together with the total source 
function are plotted in Figs.l and 2. 



A. 2. Gauge boson dominated annihilation 

In addition to the electron produced by the direct decays of vector bosons (typical 
energy of ~ there is a continuum spectrum of electrons at £■ ^ -^x/^ which are 

produced as secondary decay products {W^ ^ ^ e^, c ^ e^) and from the 

decays of charged pions produced in hadronic decays of the and bosons. Since the 
electron distributions arising from the W and the Z bosons are qualitatively similar, we 
will restrict our analysis to the W case. Following KT, the source distribution function of 
secondary-decay electrons, integrated over the quark and lepton energies, is 

Qw,^l,T,c{E, r) = n\{aV) ACw,^l,T,cfw,^i,T,c{E) (A13) 
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where 



Cw,iJ,,T,c — Bw^n + 0.17 Bw^T + 



(A14) 



and 



fw,IX,T,ci^) - ' 



—In 

M^S''"' 



(1-5) 

kM^{l+S) 
2E 



ioT E <^{1 + S), 

,'il + 5)<E<'^{l + 5), 



for 







iorE>^{l + 6) 



(A15) 



In this expression, 5 — [1 — (Mw/M^)^]^^'^ with Mw — 80 GeV being the W boson mass. 
Adopting values -Biy^M = 0.11, Bw^t — 0.11 and -Biy^c — 0.34 for the previous branching 
ratios, we obtain Cw,f^,T,c ~ 0.1729. 

The hadronization of quarks from gauge boson decays results in a shower of charged 
pions which eventually decay into electrons (tt^ ^ /x^ ^ e''^). Following KT, the source 
function for the pion-produced electrons, integrated over the quark energy distribution, 
writes 



(A16) 



where (^w,tv ~ 2/3 and 

fw,7v{E) 



M^d Je„ 



dE' 



93g-68£;/B'^5gg-27.6£/B' 



(A17) 



where Emm ^ - 5)/2 and E^ax 

case writes 

Qww{E, r) = Qw,n,T,c{E, r) + QwA^^ f) 
and is compared to the source function Qff{E,r) in Figs.l and 2. 



M^(l + 5)/2. The total source function for this 

(A18) 
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Fig. 1. — The source spectrum for a model in which fermions dominate the neutralino 
annihilation products {xx f f)'- fii'st generation, prompt electrons (PI), second generation, 
prompt electrons (P2) and secondary electrons produced in the decay of charged pions (vr). 
The light solid curve is the total source spectrum for this case. 

We also show the source spectra for the model in which gauge bosons dominate the 
annihilation. The contributions from the decays (W^ — > e^, ^ c ^ e^) and from 
the decays of charges pions produced in hadronic decays of and bosons are labelled 
KTi (light dotted curve) and KT^^ (light dashed curve), respectively. The hght solid curve 
is the total source spectrum for the case {xx ~^ WW). A neutralino density n-^ = 1 cm~^ 
and annihilation cross section {(tV)a_ = 10^^'' cm'^ s^^ have been used in this plot. 




Fig. 2. — The quantity E^Qe{E,r) is shown for the source spectra plotted in Fig.l. Curves 
are labelled as in Fig.l. 
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Fig. 3. — The overall equilibrium spectrum Ne{E,r) of the electrons whose source functions 
are shown in Fig. 1. The equilibrium spectra are calculated separately for each source 
electron population given in Fig.l. The heavy and light solid curves are the total equilibrium 
spectra for the cases XX f f and XX ~^ WW, respectively. A neutralino with mass 
My. = 100 GeV has been adopted. 
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Fig. 4. — The Coma radio halo spectrum predicted in a model in which x annihilation is 
dominated by fermions (xx ~^ f f)- The curves are for a uniform magnetic field of Buniform = 
1.3/iG (dashed) and for a radially decreasing magnetic field with central value Bq = SfiG 
(solid). We also show the radio halo spectra obtained from a power-law approximation, 
Qe ~ E~^-^ to the true source spectrum shown in Fig.l with Buniform = 1-3/iG (dot-dashed) 
and Bq = 8fiG (dotted). A constant core profile with = 0. 4/2,50^ Mpc, P = 0.76 has been 
adopted. The radio halo emission has been integrated out to l-3h^Q Mpc. Data are taken 
from Deiss et al. (1997). 
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Fig. 5. — The Coma radio halo spectrum predicted in a model in which x annihilation 
is dominated by fermions {xx ~^ f f'l solid curve) is compared to the model in which 
annihilation is dominated by higgsinos {xX ~^ WW; dotted curve). The power-law 
approximation for the case xX ~^ ff is also shown (dashed curve). The curves are plotted 
for a radially decreasing magnetic field with central value Bq = 8/iG. A constant core profile 
with Tc = OAh'^Q Mpc, /3 = 0.76 has been adopted. The radio halo emission has been 
integrated out to 1.3/i^o^ Mpc. Data are taken from Deiss et al. (1997). 
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Fig. 6. — In this figure we show how the uncertainties in the electron source spectrum reflect 
on the shape of the radio halo spectrum of Coma. Shaded areas enclose radio halo spectra 
evaluated considering both an uncertainty of a factor ±2 in normalization and ±10% in the 
slope of the source spectra. Dark gray region refers to the XX ^ ff case and pale gray 
region to the xX ^ WW case. Curves are labelled as in Fig. 5. 
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Fig. 7. — The radio halo spectrum of the cluster 1E0657-56 in the inner (empty circles) 
and outer (filled circles) regions. Dashed and dot-dashed curves are for a uniform field 
Buniform = 2yuG and Buniform = 9/iG, respectively. Solid and dotted curves are for a 



1 + {r/rc,Bf 



-0.5 



with a central amplitude 



magnetic field with a radial dependence B = Bq 
Bq = 100/iG and Bq = 90yuG, respectively. Curves are shown for a model in which neutralino 
annihilation is dominated by fermions {xx ~^ ff)- Constant core density profiles have 
been adopted here with parameters Rhaio = '^h'^o Mpc, Tc = O.SSh^Q Mpc, /5 = 0.7 and 
%,o = 9 ■ 10~^ cm~^ for the outer region and Rhaio = O-^h^o Mpc, Vc = O.OSh^Q Mpc, 
f3 = 0.49 and n^^o = 9 ■ 10^^ cm^^ for the inner region (see Liang et al. 2000). Data are from 
Liang et al. (2000). 




Fig. 8. — The Coma radio halo spectrum for a constant core beta-profile (left panel) and for 
a NFW profile (right panel). Dashed and solid curves refer to the Buniform = 1-3/i G and 
Bq = 8fiG for the beta-profile (left panel) and to Buniform = 0.45/i G and Bq = lAfiG for 
the NFW profile (right panel). A neutralino annihilation model XX ~^ ff has been assumed. 
Data are from Deiss et al. (1997). 
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Fig. 9. — The J14 — T correlation is shown for the radio halo clusters aX z > 0.2 (dark dots) 
and at z < 0.2 (gray dots). The dotted line is the best fit power-law to the data (the shaded 
area contain the uncertainty region of the best fit points, see Section 6 for details). The heavy 
lines are the relations J1.4 ~ 2"4.25 gxpected in the neutralino annihilation model discussed 
in the paper and are evaluated at z = (solid) and at z = 0.25 (dashed), respectively, data 
are taken from Feretti (1999), Giovannini et al. (1999), Liang (1999), Owen et al. (1999). 
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Fig. 10. — The gamma-ray emission predicted from XX annihilation in Coma (sohd curve) 
is compared to the sensitivities of the operating and planned gamma ray experiments. 
A neutralino with mass = 100 GeV has been assumed. The shaded area shows 
the uncertainty in the gammma-ray spectrum due to an uncertainty of a factor 2 in the 
normalization of the vr" source spectrum (see Section 3 for details). The short and long 
dashed curves refer to the gamma ray emission produced by relativistic bremmstrahlung of a 
population of primary cosmic rays and are calculated using the formula given in Sreekumar et 
al. (1996) for two choices of the Coma magnetic field B = 0.3/uG and B = l/xG, respectively. 
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Fig. 11. — The gamma-ray emission predicted from XX annihilation in galaxy clusters 
(solid curve) calculated for different neutralino masses: = 70 GeV (dot-dashed curve), 
= 100 GeV (solid curve), = 200 GeV (dotted curve) and = 300 GeV (dashed 
curve). The sensitivities of the operating and planned gamma ray experiments are also 
shown. Combined gamma ray observations of clusters from ~ 1 GeV to 500 GeV can clearly 
determine the neutralino mass from both the intensity of the spectrum and its high-energy 
cutoff (see Section 7 for details). 



